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Andreev reflection, which corresponds to the tunneling of two electrons from a metallic lead to a 
superconductor lead as a Cooper pair (or vice versa), can be exploited to measure high frequency 
noise. A detector is proposed, which consists of a normal lead-superconductor circuit, which is 
capacitively coupled to a mesoscopic circuit where noise is to be measured. We discuss two detector 
circuits: a single normal metal - superconductor tunnel junction and a normal metal separated 
from a superconductor by a quantum dot operating in the Coulomb blockade regime. A substantial 
DC current flows in the detector circuit when an appropriate photon is provided or absorbed by 
the mesoscopic circuit, which plays the role of an environment for the junction to which it couples. 
Results for the current can be cast in all cases in the form of a frequency integral of the excess noise 
of the environment weighted by a kernel which is speciflc to the transport process (quasiparticle 
tunneling, Andreev reflection,...) which is considered. We apply these ideas to the measurement 
of the excess noise of a quantum point contact and we provide numerical estimates of the detector 
current. 



I. INTRODUCTION 



Over the last decades, the measurement of noise has 
become a widely accepted diagnosis in the study of elec- 
tronic quantum transporti*^!^. Indeed, noise provides in- 
formation on the charge of the carriers in unconventional 
conductors, when considering the Fano factor - the ra- 
tio of the noise to the average current - in the regime 
where the carriers tunnel independently. Away from this 
Poissonian regime, noise also contains crucial information 
on the statistics of the charge carriers^. Experimentally, 
low frequency noise in the kHz-MHz range is more ac- 
cessible than high frequency (GHz-THz) noise: it can in 
principle be measured using state of the art time acqui- 
sition techniques. For higher frequency measurements, it 
is becoming necessary to build a noise detector on chip 
for a specific range of high frequencies. In this work we 
consider a detector circuit which is capacitively coupled 
to the mesoscopic device. This circuit is composed of a 
normal metal - superconductor junction (NS junction). 
Transport in this circuit occurs when the electrons tun- 
neling between the normal metal and the superconductor 
are either tunneling elastically, or when they are able to 
gain or to loose energy via photo-assisted Andreev re- 
flections. The "photon" is provided or absorbed from 
the mesoscopic circuit which is capacitively coupled to 
the NS detector. The measurement of a DC current in 
the detector can thus provide information on the absorp- 
tion and on the emission component of the current noise 
correlator. 

Several theoretical efforts have been made to describe 
the high frequency noise measurement process^. This is 
motivated by the fact that in specific transport setups, 
high frequency noise detection is required in order to 
fully characterize transport. Examples are Bell inequal- 
ity tests^ii and the detection of the anomalous charges 
in one dimensional correlated systems^iS. The former re- 



quire information about high frequency noise in order to 
establish a correspondence between electron number cor- 
relators and current noise correlators. The latter requires 
a high frequency noise measurement in order to obtain a 
non-zero signal when the nanotube is connected to Fermi 
liquid leads, which tend to wash out the features of a one 
dimensional correlated electron system. 

In Ref. [Tol . a detection circuit, which was capacitively 
coupled to the mesoscopic circuit to be measured, was 
proposed as a high frequency noise detector. This inter- 
esting theoretical idea has so far eluded experimental ver- 
ification, possibly because in the double dot system, ad- 
ditional (unwanted) sources of inelastic scattering render 
a precise noise measurement quite difficult. It is therefore 
necessary to look for detection circuits which are less vul- 
nerable to dissipation, as is the case of superconducting 
circuits, because of the presence of the superconducting 
gap. Indeed, the basic idea of Ref. [13 was implemented 
experimentally reccntlyiiii^ using a superconductor - in- 
sulator - superconductor (SIS) junction as a detector, 
measuring in this case the finite frequency noise charac- 
teristics of a Josephson junction. 

More recently, the noise of a carbon nan- 
otube/quantum dot was also measuredi^ using ca- 
pacitive coupling to an SIS detector, with the detection 
of Supcr-Poissonian noise resulting from inelastic cotun- 
neling processes. The latter proposals, together with 
their successful experimental implementations, indicate 
that superconducting detector circuits have advantages 
over normal detection circuits because of the presence of 
the superconducting gap. 

The purpose of the present work is thus to analyze 
a similar situation, except that the SIS junction is re- 
placed by an NS circuit which transfers two electrons 
using Andreev reflection between a normal lead and a 
superconductor. The present scheme is similar in spirit 
to the initial proposal of Ref. [l^, in the sense that it ex- 



FIG. 1: Schematic description of the set up: the meso- 
scopic device to be measured is coupled capacitively to the 
detector circuit. The latter consists of a normal metal lead- 
superconductor junction with a DC bias. 

ploits dynamical Coulomb blockade physics^. However, 
here two electrons need to be transferred from / to the 
superconductor, and such transitions involve high lying 
virtual states which are less prone to dissipation because 
of the superconducting gap, similarly to the SIS detector 
of Refs. Illlll2lll3l . Andreev reflectior>i^ typically assumes 
a good contact between a normal metal and a supercon- 
ductor, but in general it can be applied to tunneling con- 
tacts - it then involves tunneling transitions via virtual 
states. Consequently, depending on the applied DC bias, 
two successive "inelastic" electron jumps are required for 
a current to pass through the measurement circuit. The 
amplitude of the DC current as a function of bias voltage 
in the measurement circuit provides an effective readout 
of the noise power to be measured. 



A. Detector consisting of a single NS junction 

The detector circuit is depicted in Fig. [TJ Two ca- 
pacitors are placed, respectively, between each side of 
the mesoscopic device and each side of the NS tunnel 
junction. This means that a current fluctuation in the 
mesoscopic device generates, via the capacitors, a volt- 
age fluctuation across the NS junction. In turn, the volt- 
age fluctuations translate into fluctuations of the phase 
around the junction. The presence of the neighboring 
mesoscopic circuit acts as a specific electromagnetic en- 
vironment for this tunnel junction, which is described 
in the context of dynamical Coulomb blockadei^ for this 
reason. 

Fig. [2] depicts several scenarios for transport through 
an NS interface. Elastic transfer of single electrons can 
occur if the voltage applied to the junction is larger than 
the gap (Fig. [5^). Below the gap, elastic transport can 
only occur via Andreev reflectioi*i^, effectively transfer- 
ring two electrons with opposite energies with respect to 
the superconductor chemical potential (Fig. [2)d). Single 
electron can be transferred with an initial energy below 
the gap, provided that a "photon" is provided from the 
environment in order to create a quasiparticle in the su- 
perconductor (Fig. [^t). Similarly, Andreev reflection can 
be rendered inelastic by the environment: for instance. 




c) d) 




FIG. 2: Electronic tunneling in a NS junction: a) Quasipar- 
ticle electron tunneling, b) Andreev reflection, c) Photo- 
assisted electron tunneling as a quasiparticle in the supercon- 
ductor, d) Photo-assisted Andreev reflection. 



two electrons on the normal side, with total energy above 
the superconductor chemical potential, can give away a 
"photon" to the environment, so that they can be ab- 
sorbed as a Cooper pair in the superconductor (Fig. HJi). 
As we shall see, such inelastic Andreev processes are par- 
ticularly useful for noise detection. 



B. Detector consisting of an NS junction separated 
by a quantum dot 

After studying the noise detection of the single NS 
junction, we will turn later on to a double junction con- 
sisting of a normal metal lead, a quantum dot operating 
in the Coulomb blockade regime, and a superconductor 
connected to the latter (Fig. [3]). The charging energy of 
the dot is assumed to be large enough that double oc- 
cupancy is prohibited. This setup has the advantage on 
the previous proposal that additional energy flltering is 
provided by the quantum dot. Below, we refer to this 
system as the normal metal-dot-superconductor (NDS) 
detector circuit. In Fig. [H the quantum dot level is lo- 
cated above the superconductor chemical potential, and 
placed well within the gap in order to avoid quasipar- 
ticle processes. Because double occupancy is prohibited 
by the Coulomb blockade, Andreev transport occurs via 
sequential tunneling of the two electrons. Yet, because of 
energy conservation, the same energy requirements as in 
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FIG. 3: Schematic description of the NDS set up: the meso- 
scopic device to be measured is coupled capacitively to the 
detector circuit. The latter consists of a normal metal lead - 
quantum dot - superconductor junction with a DC bias. 



Fig. [2)d have to be satisfied for the final states (electrons 
with opposite energies, see Fig. [4^). In T-matrix ter- 
minology, for this transition to occur, virtual states cor- 
responding to the energy of the dot are required, which 
suppress the Andreev tunneling current because of large 
energy denominators in the transition rate. Figs. |4]b, c, 
d, describe the cases where an environment is coupled to 
the same NDS circuit. Provided that this environment 
can yield or give some of its energy to the NDS detector, 
electronic transitions via the dot can become much more 
likely because electron energies on the normal side can be 
close to that of the dot level. Such transitions can thus 
occur even if the chemical potential of the normal lead 
exceeds that of the superconductor. As we shall see later 
on, the bias voltage can act as a valve for photo-assisted 
electron transitions. It is precisely these latter situations 
which will be exploited in order to measure the noise of 
the measuring circuit (the "environment"). 

The paper is organized as follows. The model and the 
calculation of the photo-assisted current in the NS junc- 
tion are given in Sec. II, and results for the detection of 
the noise of a quantum point contact arc presented. In 
Sec. Ill, we perform the same analysis for the NDS junc- 
tion and compare the results. The conclusion is point 
out in Sec. IV. 



II. TUNNELING CURRENT THROUGH THE 
NS JUNCTION 



A. Model Hamiltonian 

The Hamiltonian which describes the decoupled 
normal metal lead-superconductor-environment (meso- 
scopic circuit) system reads 
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FIG. 4: Andreev reflection in the NDS junction: a) Andreev 
reflection in the elastic regime, b), c) Photo-assisted Andreev 
reflection, where a "photon" is provided to or provided by 
the environment, d) Absorption of a Cooper pair with (re- 
verse) photo-assisted Andreev reflection, where a "photon" 
is provided by the environment. For cases b), c), d), which 
require passing through the dot, the tunneling of electrons is 
sequential. 



where 



(2) 



k,{7 



describes the energy states in the lead, with ^ an elec- 
tron creation operator. The superconductor Hamiltonian 
has the diagonal form 



Hos - f^sNs = J2 ^1^1 



(3) 



where 75,0-, 7^ a are quasiparticle operators, which relate 
to the Fermi operators c^.o-, c^^- by the Bogoliubov trans- 



formation 



C-q.l = Uql-q,l " "97g,t ' 



(4) 



and E„ 



+ Cn is the quasiparticle energy, Q 



Ho — Hq + Hqs 



Hp 



(1) 



Eg — /^s is the normal state single-electron energy counted 
from the Fermi level /xg, A is the superconducting gap 
which will be assumed to be the largest energy scale in 
these calculations. Hereafter, we also define eV = hl — 
fis a-nd assume fis = 0. 
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Here we do not specify the Hamiltonian of the envi- 
ronment because the environment represents an open sys- 
tem: the mesoscopic circuit which represents the environ- 
ment will only manifest itself via the phase fluctuations 
{(j){Q)(t){t)) which are induced at the NS junction (or later 
on for the NDS circuit, at the dot-superconductor junc- 
tion) because of the location of the capacitor plates. In 
what follows, wc shall assume that the unsymmetrized 
noise spectral density 

/ + 00 
dte*"*((/(0)/(t))) , (5) 
-oo 

corresponding to photon emission (for positive fre- 
quency), or alternatively iS'~(a;) = S~^{—lu), the spectral 
density of noise corresponding to photon absorption, arc 
both specified by the transport properties of the meso- 
scopic circuit^ i^^'^^ . Here ((. . .)) stands for an irreducible 
noise correlator, where the product of average currents 
has been subtracted out. 

The tunneling Hamiltonian describing the electron 
transferring between the superconductor and normal 
metal lead in the NS junction is 

Ht=Y1 TkA^aV-'" > (6) 

the indices fc, q refer to the normal metal lead, supercon- 
ductor. Wc consider for simplicity that Tk^q = Tq. Note 
that the tunneling Hamiltonian contains a fluctuating 
phase factor, which represents the coupling to the meso- 
scopic circuit. Indeed, because of the capacitive coupling 
between the sides of the NS junction and the mesoscopic 
circuit, a current fluctuation translates into a voltage 
fluctuation across the NS junction. Both are related via 
the trans-impedance of the circuiti^ ^(^) Z{uj)I{u!). 
Next, the voltage fluctuations translate into phase fluctu- 
ations across the junction, as the phase is the canonical 
conjugate of the charge at the junctionsii: the phase 
is thus considered as a quantum mechanical operator 
throughout this paper. For definition purposes, it is con- 
venient to introduce the correlation of the phase opera- 
tors 

jit) = mt)-m]m) , (7) 

where the phase operator is related to the voltage by the 
relation 

0(t) = e I dt'V{t') . (8) 

Given a specific circuit (capacitors, resistances,...) the 
phase correlator is therefore expressed in terms of the 
trans-impedance of the circuit and the spectral density 
of noisei^ 

J(t) = |^pa;^5,(c.)(e---1), (9) 



where Rk = 27r/e^ is the quantum of resistance and 
Si[oj) = S+{-ijj). 

The present system bears similarities with the study of 
inelastic Andreev reflection in the case where the super- 
conductor contains phase fluctuation a^^'^^ . Such phase 
fluctuations destroy the symmetry between electrons and 
holes, and affect the current voltage characteristics of the 
NS junction. 

B. Tunneling current 

The tunneling current associated with two electrons is 
given by the Fermi Golden Rule / = 2eTi^f, with the 
tunneling rate^^ 

r,^/=2^^|(/|T|z)|^<5(6.-e;) , (10) 
/ 

where and ef are the tunneling energies of the ini- 
tal and final states, including the environment. T is the 
transition operator, which is expressed as 

T = Ht + HtY.(-. tt—Ht] , (11) 

n—l ^ ^ 

with ?7 is a positive infinitesimal. 

Throughout this paper, one considers the photo- 
assisted tunneling (PAT) current due to the high fre- 
quency current fluctuations of the mesoscopic device, as 
the differenced^: 

IpAT = /(environment) — /(no environment) , (12) 

where, in general, the total current for tunneling of elec- 
trons through the junction is / = /^ — /<_. 

However, experimentallyi^ it is difficult to couple ca- 
pacitively and then to remove the mesoscopic device cir- 
cuit from the detector circuit. What is in fact often mea- 
sured is the excess noise, i.e. the difference between cur- 
rent fluctuations at a given bias and zero bias in the 
mesoscopic circuit. In the latter on, we will calculate 
the excess noise S'+,(w) of S^{uj). In this work, we thus 
measure the difference between the currents through the 
detector when the device circuit is applied a bias volt- 
age and zero bias, as a function of detector bias voltage, 
which is defined as 

AIpATieV) = lieVd + 0, eV) ~ I[eVd = 0, eV) . (13) 

This also corresponds to the difference between the PAT 
currents through the detector when the device circuit 
is applied a bias voltage and zero bias A/p^T(el/) ~ 
IpAxieVd 7^ 0,eV) — IpAxieVd = 0,6^) because the 
contributions with no environment cancel out. The dif- 
ference PAT current thus provides crucial information on 
the spectral density of excess noise of the mesoscopic de- 
vice. Notice that our calculation applies to the zero tem- 
perature case for convenience, but it can be generalized 
to finite temperatures. 
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C. Single electron tunneling 

Although our focus of interest concerns photo-assisted 
Andreev reflection, we need to compute all possible con- 
tributions. The current associated with one electron tun- 
neling is given by the Fermi golden rule: 



/ = 2^e^|(/|i/TN)|'<5(e.-e/) 
/ 



(14) 



The calculation of the current proceeds in the same 
way as that of a normal metal junctioni^, except that one 
has to take into account the superconducting density of 
states on the right side of the junction which is done by 
exploiting the Bogoliubov transformation. For the case 
of electrons tunneling from the normal metal lead to the 
superconductor {eV > A), the current from left to right 
reads 



' -oo J A 



Rk 



de 

V 



dE- 



E 



VE^ - A2 
E 



27r 

|2 



dJ-^S^i.) 



S{E-e) 



de / dE , , 

A _ (e-E) 



(15) 



with Rk the quantum of resistance. This current in- 
cludes both an clastic and an inelastic contribution, the 
former being rcnormalizcd by the presence of the envi- 
ronment. Here, e is the energy of an electron in the 
normal metal lead and E is the energy of a quasipar- 
ticlc in the Superconductor lead. Changing variables in 
the inelastic term to D, = e — E, S = e + E, an using 
/ dx{x + a)/y^{x + a)^ -b^ = {x + a)^ - /2, after 
computing the current from right to left in a similar way, 
we obtain 



^IpAT = —Ci 

' — oo 

-OO 



stA-mi7\^,^v) ,(16) 



where the transmission coefficient of the NS junction in 
the normal state is defined as T = Aii^ML-N'ogT^ ^ Cie = 
eT /Rk- The weight functions arc defined as 



(17) 



Kl"/^ {n, eV) = v/(0-ey)2- A2 . (18) 

Similarly, we obtain the formula of Alp at for the case 
3V < -A 

AIpAT = -Cie ( r dJ-^^StA-^)] Kt{eV) 



|Z(-17)|2 
172 



+Cie / d^! '"-^;-" sUmiT\^.eV) . (19) 

leV+A 



For the case — A < eV < A: there are no elastic transi- 
tions because electrons cannot enter the superconducting 



gap. Nevertheless, due to the presence of the environ- 
ment, electrons can absorb/emit energy from/to the en- 
vironment so that an inelastic quasiparticle current flows. 



AL 



PAT 



/eV-A \7(Cl\\'^ 
dn^-^^stA~mT/\^,eV) 

dn'-^^^^SUmiT'i^.eV) . (20) 

eV+A " 



D. Two electrons tunneling as two quasiparticles 

The single electron (quasiparticle) tunneling current is 
of first order in the tunneling amplitude. We now turn 
to processes which invoke the tunneling of two electrons 
through the NS interface. Indeed, because our aim is 
to show that Andreev reflection can be used to measure 
noise, we need to examine all two electron processes, we 
start with the transfer of two electrons as quasiparticles 
above the gap. Calculations of the matrix element in Eq. 
(fTU|) are then carried out to second order in the tunneling 
Hamiltonian using the T-matrix: 



I = A-kcY^MHt 
f 



1 



hj - Ho + e, 



■HT\t)\ S{e, -ef) , (21) 



The initial state is a product: 



\i) = \Gl)<E>\Gs)<E>\R) 



(22) 



where |Gl) denotes a ground state, which corresponds to 
a filled Fermi sea for the normal electrode. \Gs) is the 
BCS ground state in superconductor lead. \R) denotes 
the initial state of the environment. On the other hand, 
our guess for the final state should read 



I/) - AA',^'<-i',^'\^^) ® ® \R') 



(23) 
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when two electrons are emitted from superconductor. 
The "guess" of Eq. (|23p is an informed one: an elec- 
tron pair is broken in the superconductor, and one elec- 
tron tunnels to the normal metal lead, while the other 
becomes a quasiparticle in the superconductor; the same 
process for the second electron tunneled to the normal 
metal lead. When superconductor lead absorbs two elec- 
trons 

|./> = Ck,.cu:,,nU^l,^^, \Gl) ® \Gs) ® , (24) 

The "guess" of Eq. and ((24|) is: two electrons can tunnel 
from the normal metal lead and become quasiparticles in 
the superconductor. Here, \R') is the final state of the 
environment. 

Introducing the closure relation for the eigenstates of 
the non-connected system and using the fact that 

{v\{ei-HQ±iri)~'^\v) = =Fi dte*^'''"''"^"')*, one can ex- 
ponentiate all the energy denominators. Then, by trans- 
forming the time dependent phases into a time depen- 
dence of the tunneling Hamiltonian, we can integrate 
out all final and virtual states. This allows to rewrite 
the tunneling current in terms of tunneling operators in 
the interaction representation, to lowest order 0(Tq). 

For the case of two electrons tunneling from Supercon- 
ductor lead to normal metal lead, the current reads 

/OO /"OO /"OO 

dt / dt' / dt"e-"(*'+*")e*(^^-*'^)(2*-*'-*") 
-OO JO Jo 

X (i?^(i - t")Hl{t)HT{t')HT{0)) . (25) 
Here we consider quasiparticle tunneling so 

(cUi - t"H2a, (Ocg3a3 (i')c?4^4 (0)) 

+ |«,n„^fe-*^«(*-*")e-'^.'(*-*') . (26) 



The exponentials of phase operators are calculated with 
the definition of Eq. so tha,^^ 



^gJ</.(t-t")g#(t)g-i</'(i')g-#<) 
_ J{t~t"-t') + J{t-t") + J{t~t') + J{t)~J{~t")-J{t') 



(27) 



then we further assume that J{t) <^ 1, which means a low 
trans-impedance approximation together with the fact 
that J{t) is well behaved at large times. This allows 
to expand the exponential of phase correlators e'''-*'' « 
1 -I- J{t). The result for the current contains both an 
elastic and an inelastic contribution = + T^'^'- ^ 
where 

(28) 

with "fo^, ^'1^^(01, eV, rj) are defined as in Eqs. ([55)1 and 
([54)) in Appendix A, and 

dn^ItS:}lsj{-{n))iq:t{n,eV,ii) , 



Wtt'^Rk J2A+2. 



2eV 



f72 



where K2^^{n, eV, rf) is defined in Eq. ([55)) in Appendix 
A. The elastic contribution exists only if eV < —A. For 
—A < eV only the inelastic part contributes to I^. 

Similarly we calculate for J_>. There is a symmetry 
between the magnitude between the right and left moving 
current upon bias reversal: the expression for /<_ is the 
same as I—,, if we replace —eV by eV . 

So in the interval \eV\ < A, we obtain 



AIpATieV)= C2e / dn 



l2A-2eV 

-C2e r dn^-^tI}}ls+^{n)KrS . (30) 

J2A+2eV ""^ 



with K'2"et{eV) 



Kllti-aV) and C2e = eT^/lGn^RR. 



E. Two electron tunneling as a Cooper pair: 
Andreev reflection 



In this case, we also needs to carry out calculations of 
the matrix element in Eq. ([TO)) to second order in the 
tunneling Hamiltonian. Typically, the initial state will 
be as shown in Eq. ([22)) . On the other hand our guess 
for the final state reads: 



cI_„\\Gl)®\Gs)®\R') , 



(31) 

when a Cooper pair is emitted from superconductor, or 

I/) = 2-'/^[ck.aCk'.-a - Ck',^Ck,-a]\GL) ® \Gs) ® \R') , 

(32) 

when Superconductor lead absorbs a Cooper pair. Here, 
\R') is the final state of the environment. The "guess" 
of Eqs. ([5T)) and ([5^ is again an informed one: indeed, 
the s-wave symmetry of the superconductor imposes that 
only singlet pairs of electrons can be emitted or absorbed. 
This phenomenon has been described in the early work 
on entanglement in mesoscopic physics^ii^, and the re- 
sulting final state can in principle be detected through a 
violation of Bell inequalities^,. For this Andreev process, 
we can now write the tunneling current as a function 
of the normal (and anomalous) Green's functions of the 
normal metal lead, Glo-, the quantum dot, Gdu, and the 
superconductor, F„ (see Appendix B), which is the same 
as in Ref. [H 



dt 



dt' / dt"e-'K*'+*")e'(Ms-M.)(2t-t'-t") 



X E {-GUk.t-t" -t')G>^„{k\t)F;{<^^^t")F^,{q,t') 



E 

k.k' ,q.q' .(7 

+ GL(fc, t - t")G>_,(fc', t ~ t')F:{q', ^t")F,{q, t')] 

^^J(t~t"-t') + J(t-t") + J(t-t') + J{t)-J(-t")-J(t') _ j-gg-j 

The result for the current contains both an elastic and 
an inelastic contribution: 



(29) 



Fl + I: 



inel 



(34) 
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where the elastic contribution reads 



e72 
27r3 



de I dE 



eV 



dE' ^ — 

A y/E^ - A^VE'^ - A2 



1 - 



47r 



1 



27r 



The inelastic contribution to is 

/T-2 /'OO pOC pOC pQC 

~ / de / de' / dE / di;' 

-flK JeV JeV J A J A 

A2 |^(_(e + e'))p5,(-(e + e')), 



(35) 



V£^2 _ A2^^/2 (e + e')2 



(36) 



where the denominators are specified in Appendix C. /_> 
is derived in a similar manner, but its expression is omit- 
ted here. Nevertheless, it effects will be displayed in the 
measurement of the noise of a point contact. 

The above expressions constitute the central result of 
this work: one understands now how the current fluctu- 
ations in the neighboring mesoscopic circuit give rise to 
inelastic and elastic contributions in the current between 
a normal metal and a superconductor. 

We find that both current contributions have the same 
form and the current fluctuations of the mesoscopic de- 
vice affect the current in the detector at the energy cor- 
responding to the total energy of two electrons in the 



J 



normal lead. Andreev reflection therefore acts as an en- 
ergy filter. 

Next, we can change variables as in the previous sec- 
tions. With an arbitrary bias eV, we obtain 



AIpATieV) 

= —Cms 
Cns 



+00 \7(uj)\^ 

dLu^-^J^St....i~u;)K'Nsi^, el/, rj) 



|2 



dVL'- 



c 



NS 



rf» '^L?'' 5+ ,,,,(»)gj,"g'(f}, eV, 77) , (37) 



2eV 



with Cns = eT^A'^f-K^Rx- The kernel fimctions Kf^g, 
K'^n's shown in Appendix C. 



F. Quantum Point Contact as a source of noise 

1. Spectral density of excess noise 

In this section, we illustrate the present results with 
a simple example. We consider for this purpose a quan- 
tum point contact, which is a device whose noise spec- 
tral density is well characterized by using the scattering 
theoryii2^. Here however, we consider unsymmetrized 
noise correlators: 



5+(a>) 



2^T(1 - T){eVd - nu)e{eVd - htu) , 
'-[-2T^nuj - r(l - T){eVd + hLu)0{-eVd - huj) + T{1 

I 



if w > 0, 

T){eVd~huj)] , if < , 



(38) 



where is the Heaviside function and T is the trans- 
mission probability, which is assumed to have a weak 
dependence on energy over the voltage range eV . 

As pointed out above, we are computing the difference 
of the PAT currents in the presence and in the absence 
of the DC bias. This means that we insert the spectral 
density of excess noise of the mesoscopic device, which 
for a point contact bears most of its weight near zero 
frequencies. Excess noise decreases linearly to zero over 
a range [0, ±eVd] for positive and negative frequencies: 

StA^)^—T{l - T){eVd - \nuj\)e{eVd - \hLo\) , (39) 



2. Numerical calculations 

We choose a generic form for the transimpcdancc, sim- 
ilar in spirit to that chosen in Ref. [lol . Considering the 



circuit in Figs. [T]and[3l aX uj ~ Q, device and detector 
are not coupled and the transimpedence should therefore 
vanish. On the other hand, the transimpedance is pre- 
dicted to have a constant behavior at large frequency. 
We therefore choose the following generic form for the 
transimpedance 



(40) 



where R is the typical high frequency impedance and the 
crossover frequency is estimated from the experimen- 
tal data of Ref. [12, choosing a finite cutoff frequency ujq 
means that at frequencies w <C wq, the mesoscopic cir- 
cuit has no influence on the detector circuit because low 
frequencies do not propagate through a capacitor. 

We calculate numerically the PAT currents in the three 
above cases: single and two quasiparticle tunneling, and 
Andreev reflection. All energies are measured with re- 
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FIG. 5: A7pAT plotted as a function of the DC bias voltage, for a mesoscopic device voltage bias eVd- 0.3 (continuous line), 
0.5 (dashed line) and 0.8 (dotted hne). The left /center /right panels depict single quasiparticle tunneling/two quasiparticle 
tunneling/ Andreev reflection. A/pat is in units of C2e, with T = 0.6 (see text). 



spect to the superconducting gap A. In these units we 
chose ujQ = 0.3 and 77 = 0.001. Currents are typically 
plotted as a function of the DC bias voltage eV of the 
detector, for several values of the mesoscopic bias voltage 
eVd (the "environment" ) . Our motivation is to consider 
the PAT currents with the condition \eV\ < 1 {\eV\ < A), 
where the effect of the environment on the PAT current 
is most pronounced, and we shall predict that two elec- 
trons tunneling as a Cooper pair (Andreev processes) 
contributes the most to the PAT current, except close 
to eV = A. Because of the symmetry between negative 
eV and positive eV, we display the results for eV > 0. 
The PAT currents for the three above processes are plot- 
ted next to one other in Fig. [5] for comparison. 

We find in Fig. [5] that if the chemical potential of the 
normal metal lead is close to the potential of the super- 
conductor lead {eV ^ A), the PAT currents in the three 
cases are suppressed. For the two cases of quasiparti- 
cle tunneling, the PAT currents are equal to zero below a 
certain threshold. The single quasiparticle current differs 
from zero at a threshold which is identified as A — eVd, 
that is quasiparticles are able to tunnel above that su- 
perconductor gap only if they can borrow the necessary 
energy from the mesoscopic device. This explains why 
the curves associated with different values of the meso- 
scopic device bias voltage are shifted to the right as eVd 
decreases. For two quasiparticle particle tunneling, we 
observe that the PAT current has a similar threshold, 
which, compared to Fig. [5^, is pushed toward the right 
in Fig. ]Sjp, because more energy is needed to transfer two 
electrons above the gap, compared to a single electron. 
Nor surprisingly, the curves corresponding are once again 
shifted to the right with decreasing eVd- These curves all 



have a sharp maximum at eV ~ A. 

We turn now to the Andreev PAT difference current, 
which dominates the two above processes at small and 
moderate bias. Note that the total Andreev current con- 
tains an elastic contribution as well as an inelastic con- 
tribution below the gap, contrary to quasiparticle tun- 
neling which have contributions below the gap only be- 
cause these processes are photo-assisted. Because we are 
computing the difference between PAT current with and 
without the mesoscopic bias voltage, we expect that the 
clastic contribution cancels out. However the first term 
of Eq. ([37|) tells us that the presence of the environ- 
ment also gives rise to an effective elastic contribution 
to the PAT Andreev difference current. Unfortunately, 
this elastic correction is not small compared to the true 
Andreev current. 

Looking at Fig. [SJ we note that the PAT curve for 
Andreev processes is shifted to the left when the bias of 
the mesoscopic device is increased. The environment pro- 
vides/absorbs energy to pairs of electrons whose energies 
are not symmetrical with respect to the superconducting 
chemical potential. At very weak eV, the elastic cor- 
rection is small compared to the true Andreev current. 
We expect the PAT current contributions to originate 
from pairs of electrons in the normal metal below the 
superconducting chemical potential which can extract a 
photon from the environment. At the same time. Cooper 
pairs can be ejected in the normal metal as a reverse An- 
dreev processes provided they borrow a photon to the 
environment. There is a balance between the right and 
left currents at very weak bias. 

As the bias is increased, electrons pairs whose ener- 
gies are above this chemical potential will now be able to 
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yield a photon to the environment, giving rise to an in- 
crease of the inelastic PAT current. Also, the reverse An- 
dreev process mentioned above becomes more restricted 
because the available empty states for electrons in the 
leads lie higher at positive bias. The clastic PAT current 
(not shown in the figures) increases when we increase the 
detector bias but it is always smaller than the inelastic 
contribution. The total PAT current increases gradually 
(Fig. [51 right panel). 

A zoom of this Andreev contribution is made in the 
region of small bias. There is in fact no threshold for 
the PAT Andreev current: for small bias, it has a lin- 
ear behavior (Fig. (H). According to Eq. (|37p . when 
eV > eVd/2, the inelastic difference current of a Cooper 
pair tunneling from S to L by absorbing energy from 
mesoscopic device vanishes. The PAT current now de- 
scribes two electron tunneling from L to S elastically or 
inelastically. In this intermediate bias regime, the elastic 
and inelastic contributions have now a tendency to can- 
cel each other. The current reaches a maximum close to 
the gap, and then it decreases dramatically at the gap. 
This is consistent with the fact that for positive bias, the 
initial state for two quasiparticle tunnel processes and 
for Andreev reflection is precisely the same: close to the 
gap two quasiparticle tunneling takes over the Andreev 
process. It becomes more efficient for electrons to be ac- 
tivated above the gap than to be converted into a Cooper 
pair because the energy loss needed for the latter is quite 
large. Unlike a conventional normal metal - supercon- 
ductor junction with elastic scattering only, where the 
relative importance of quasiparticle tunneling and An- 
dreev reflection are interchanged precisely at the gap, 
here the dominance of quasiparticle tunneling manifests 
itself before the voltage bias reaches the gap. Note also 
in Fig. [51 that the magnitude of the Andreev current be- 
fore the two quasiparticle threshold is precisely the same 



as the magnitude of the two quasiparticle at eV = A, 
which confirms this conversion scenario. A comparison 
of the two processes is displayed in Fig. [71 

In practice, the different contributions to the PAT cur- 
rent cannot be separated: one measures the sum of the 
three contributions which arc plotted in Fig. [5l However, 
we claim that for a broad voltage range (from eV = to 
the two quasiparticle threshold), the main contribution 
to the current comes from photo-assisted Andreev pro- 
cesses. The confrontation of Eq. (|37p with an experimen- 
tal measurement of the PAT current below the gap could 
thus serve an effective noise measurement, as the weight 
functions Kfjg{Lu,eV,T]) and K^]^g^{n,eV,ri) are known. 

Notice that in all our numerical results, the PAT cur- 
rents are plotted in units oie^R'^T'^T{l-T)A/8TT^h'^RK- 
We put some tentative numbers in these quantities. Here 
T = 0.6 is the effective transmission coefficient of the NS 
interface, A = 240/iey, T ~ 0.5 is the transmission of the 
quantum point contact to be measured, R = 0.03i?x {Rk 
is the resistance quantum) is the resistance which enters 
the transimpedance. This implies, e.g. for the PAT cur- 
rent in Fig. [6] at the top of the peak, AIpat — 10~^°A 
with the device bias Vd ~ 0.8A/e = 48/iy, which seems 
an acceptable value which is compatible with current 
measurement techniques. 



III. TUNNELING CURRENT THROUGH A 
NDS JUNCTION 

We now turn to a different setup for noise detection 
where electrons in a normal metal lead transit through 
a quantum dot in the Coulomb blockade regime before 
going into the superconductor. The essential ingredients 
are the same as the previous section, except that addi- 
tional energy filtering occurs because of the dot. In this 
section, we choose the parameters of the device so that 
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only Andreev processes are relevant. 

A. Model Hamiltonian 

The Hamiltonian which describes the decoupled nor- 
mal metal lead - dot - superconductor - environment 
(mesoscopic circuit) system reads 

Ho = Hq^ + + -ffos + Henv , (41) 

where the Hamiltonian of normal metal lead and super- 
conductor lead arc described as above. 

The Hamiltonian for the quantum dot reads 

Hod = ^'^DC^D,a^D,a + Unini , (42) 

(7 

where U will be assumed to be infinite, assuming a small 
capacitance of the dot. We consider that the dot pos- 
sesses only a single energy level for simplicity. 

The tunneling Hamiltonian includes electron tunneling 
between the superconductor and the dot, as well as the 
tunneling between the dot and the normal metal lead 

Ht = {Hti + Ht2) + h.c , (43) 

q.a 

the indices k, D, q refer to the normal metal lead, quan- 
tum dot. superconductor. Wc consider the simple case 
To.q = Ti, and Tk^o T2. 

For photo-assisted Andreev processes, one needs to 
carry out calculations of the matrix element in Eq. (jlOp 
to fourth order in the tunneling Hamiltonian. In what 
follows, we assume that the dot is initially empty, ow- 
ing to the asymmetry between the two tunnel barriers. 
The barrier between the normal metal lead and the dot 
is supposed to be opaque compared to that between the 
dot and the superconductor. As a result, the rate of es- 
cape of electrons from the dot to the superconductor is 
substantially larger than the tunneling rate of electrons 
from the normal lead to the dot (see below for actual 
numbers). 



There are two possibilities for charge transfer pro- 
cesses: a Cooper pair in the superconductor is transmit- 
ted to the normal lead or vice versa. The first process 
involves the electron from a Cooper pair tunneling onto 
the dot, next this electron escapes in the lead; the other 
electron from the same Cooper pair then undergoes the 
same two tunneling processes. Similar transitions, in the 
opposite direction, are necessary for two electrons from 
the normal lead to end up as a Cooper pair in the super- 
conductor. Note that this description of events assumes 
implicity that superconductor lead remains in the ground 
state in the initial and final states (Andreev process). On 
the other hand, if normal metal lead is initially in the 
ground state (filled Fermi sea), it is left in an excited 
state with two electrons having energies above Fermi en- 
ergy Ep in the final state. The extra energy has been 
provided by the environment. Typically, 

\i) = \Gl) ® \Gs) ® \Oqd) ® \R) , (44) 

where \G,,,) denotes a ground state, which corresponds to 
a filled Fermi sea for the normal electrode. \Oqd) is the 
vacuum of the quantum dot, and \R) denotes the initial 
state of the environment. On the other hand our guess 
for the final state should read 

I/) - ^-'^M.A'.-.-4'.A,-JGLmGsmoQD)m') , 

(45) 

when a Cooper pair is emitted from S, or 

I/) = 2-^^^CM',-a-Ck',aCk,-a]\GLmGs)(E>\0QD)®\R') , 

(46) 

when superconductor lead absorbs a Cooper pair. Here, 
\R') is the final state of the environment. The justifica- 
tion for the choice of Eqs. ([45]) and (|46|) is the same as 
in Sec. [ITeI 



B. General formula for the photo-assisted Andreev 
current 

For the case of two electrons tunneling from supercon- 
ductor lead to normal metal lead, the current reads 



/OC pOD i^OQ pOD pOD pOD pOQ 

dt / dh / dt2 / dt3 / dt[ / dt'^ / dt^e-''(*i+*^+*^+*'i+*^+*3) 
-00 Jo Jo Jo Jo Jo Jo 



X g'Ms(2i — t'l — — 2*3 — ii — 42) g — V-L (Zt — tj— 43 — 2*1 — (2 — is) 



H^^it - t[ -t'^- t'^)Hl^{t -t'^- t'.^)Hl^{t - t'^)Hl^{t)HT2{ti +t2+ h)HTi{h + t2)HT2{ti)HTim . (47) 

I 



The problem is thus reduced to the calculation of cor- Using Wick's theorem, these can be expressed in terms 
relators of the tunneling Hamiltonian in the ground state. 
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of single particle Green's function because the Hamil- 
tonian of the isolated components is quadratic (except, 
maybe for the environment, which is dealt separately). 
We can now write the tunneling current as a function 



of the normal (and anomalous) Green's functions of the 
normal metal lead, Glc^ the quantum dot, Gdu, and the 
superconductor, F^r (which are shown in Appendix B) 



/•OO /"OO /"OO /"OO 

2eT'lT^ I dt / dh / dt2 / dts / dt[ 

'-00 Jo Jo Jo Jo 



Jo 



^ gi/is(2i — t'l — — 2*3 — ii — 42) g — i/Ji (2t — — 43 — 2*1 — t2 — is) 

X [-Kii', -t'l - t'^)F-Aq, ti + ^2)GL(fc, t-t'^-t'.,-h-t2- h)G>_^{k\ t - h) 

X GL(-i'i)GL.(-i3)Gk(^3)G^-.(ti) + ^^;(9',-ti -t'2)^^.('?,ii +i2) 

X G>,(fc, t-t'^-t'^- h)G>^Jk', t-h-h- i3)GL(-^'i)Gi,_,(-t^,)G^_,(t3)G^,(ti)" 

Xe'^(*~*'l~*2-*3) + -'(*-*3) + '^(*-*i-*2-*3-*l-*2) + ^(*-t3-tl-t2)-J(-ii-t2)-'^(*l+t2) 



(48) 



The result for the current contains both an elastic and 
an inelastic contribution. The elastic contribution reads 
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(49) 



with 1)2- is the original denominator which is not af- 
fected by the environment, which is defined by Eq. ([S^ . 
and is the denominator product affected by the en- 
vironment (see Eq. (|63|1 of Appendix D). The inelastic 
contribution reads 

^,2^;2 /"OO />oo j*oo poo 

9n I de de' dE dE' 

Rk JeV JeV J A J A 

\z{-{e + e'))\'Sji-ie + e')) 



VE^ - AWe'^ ~ (e + e')^ 



where D™'^' is the denominator product attributed to 
the inelastic current, which is defined in Eq. ([64]) of Ap- 
pendix D. 

Here, 71 = 2itMsT1, 72 = 2hMnT2 define the tunnel- 
ing rates between the superconductor and the dot, be- 
tween the dot and the normal metal lead, respectively, 
with Ms and Mn the density of states per spin of the 
two metals in the normal state, at the chemical poten- 
tials /is and /iL, respectively. All contributions to the 
current contain denominators where the infinitesimal ?/ 
(adiabatic switching parameter) is included in order to 
avoid divergences. In fact, it has been shown in Refs. 
[21II26I that a proper resummation of the perturbation se- 
ries, including all round trips from the dot to the normal 
leads, leads to a broadening of the dot level. We take 



into account this broadening by replacing 77 by 7/2 with 
7 = 7i + 72 into our calculations (including only a broad- 
ening due to the superconductor). As mentioned above, 
we have assumed that 71 2> 72- In order to avoid the ex- 
citation of quasiparticles above the gap, these rates need 
also to fulfill the condition ££> -h 7 < A. In what fol- 
lows, we keep the notation rj in our expressions, bearing 
in mind that it represents the line width associated with 
the leads. For numerical purposes, it will be sufficient to 
assume that 77 is kept very small compared to the super- 
conducting gap, as well as all the relevant level energies 
(dot level position, bias voltages, etc.). 

The above expressions constitute the second main re- 
sult of this work: one understands now how the current 
fluctuations in the neighboring mesoscopic circuit give 
rise to inelastic and elastic contributions in the current 
in the NDS device. 

We find that both right and left current contributions 
have the same form. The current - current fluctuations of 
the mesoscopic device affect the detector current at the 
energy corresponding to the total energy of two electrons 
exiting in (or entering from) the normal lead. Therefore, 
we proceed to the same change of variables a for the single 
NS junction. 

For eV > 0, elastic current contributions in /_» are 
present but the same contributions in vanish (the 
opposite is true for the case of eV^ < 0). 

Changing variables in the inelastic contributions, 
and deflning the kernel functions iir|^^g(u;, eV, ed, 77), 
if ™|,'g(17, ey, e^i, 77) as in Appendix D, for eV > and 
eV < 0, we obtain: 



AIpATieV) 

= -Gnds I dij )"^^^^ S^^{-uj)K'}jjjs{uJ,eV,eD,v) 
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Jjn\^^^StAmNU^. eV, en, V) , (51) 

with Cnds = nhl^^ I'^'^Rk- 

The first term in Eq. (|5ip describes the elastic contri- 
bution in the PAT current. Although we are less inter- 
ested in this contribution, we cannot ignore it in practice 
because it contributes to the total A/p^j-. The envi- 
ronment affects this current contribution but at the end 
of the tunneling processes, there is no energy exchanged 
between the device and the detector circuit. The second 
term in Eq. (|5ip describes the tunneling of a Cooper 
pair from the normal lead to the superconductor via the 
quantum dot, with energy exchange. The electrons can 
absorb energy (in the case their total energy is smaller 
than the superconductor chemical potential or emit 
energy (if their total energy is bigger than ^s)- The last 
term in Eq. (jSip describes the inverse tunneling process: 
a Cooper pair absorb energy from the neighboring device, 
its constituent electrons then tunnel to the normal lead. 
In this event, the total energy of the outgoing electrons is 
positive. If on the contrary, this total energy is negative, 
the Cooper pair has emitted energy to the device. 

In order to understand how the detector circuit af- 
fects the behavior of the current (in the presence of 
the environment), we investigate the weight functions 
K^jjg{uj, eV,ej),ri) and K^^'j^g{^l, eV,eD, v) separately. 

The weight function Kfjjjg{u!,eV,eD,'ri) is plotted in 
Fig. [8] as a function of frequency uj, for two values of the 
bias voltage and two values of the dot level position. This 
elastic kernel is symmetric under bias voltage reversal 
{KNDsi-<^^) = -^NDsi^^'^))- From the right panel of 
this figure, where we consider \eV\ > ejy, we find that 
there is a small step at w = A — eu and a sharp peak at 
u! = —A + £/). The peak is assymetric, and its height is 
much larger than that of the step. When cu < — A + en, 
^NDS changes sign and becomes negative. The voltage 
bias eV mainly affects the amplitude of the peak and of 
the step in Kf^jjg. The left panel of Fig. [8] describes 
^NDS when \eV\ < en- The peak height decreases quite 
fast as a function of eV , and its location is shifted at 
Lj ~ —A + eV. The peak is symmetric for large bias. 

We turn now to the truly photo-assisted processes, 
which involve either absorption or emission of energy. 
The kernel K]^fjg{Q,eV,eD,'r]) is plotted in Fig. [9] as a 
function of frequency ft, which corresponds to the total 
energy of two electrons, for en > 0. In the left panel, 
eV is negative, and in the center panel, eV is positive 
but eV < ejj, and finally the right panel of Fig. [9] de- 
scribes eV > en- We find that when eV < en, there is 
a step at ^ eo + eV . When we increase eV close to 
££), the step still dominates K'^^^g but there is a small 
peak aX Vl = —A -I- eV . When eV > cd this (inverted) 
peak is very sharp. This is explicit in the right panel. 
The inverted peak, which has a large amplitude, makes 
it now difficult to observe the step. The (inverted) peak 



is located at J7 = — A 4- e^i. Again, eV mostly affects 
the ampfitude of K'^^ds- ^or < (not shown), the 
result is similar to that of e/3 > with opposite eV , but 
the amplitude of the peak is doubled compared to that 
of ££) > when \eV\ > \eD\- 

Note that understanding the behavior of the two 
weight functions as a function of the different parame- 
ters {eV and eo) of the detector is crucial. It allows to 
control the effect of the device voltage bias eVd on the 
DC current of the detector, and it is therefore the key 
for extracting the noise from the measurement of this 
DC current. 



C. Application to a Quantum Point Contact 

We now calculate AIpat from Eq. (|?T|) with the spec- 
tral density of excess noise of a Quantum Point Con- 
tact, given by Eq. ([39|) . We consider the PAT current 
as a function of the detector voltage eV for several val- 
ues of the device voltage eVd, which are shown in Fig. 
[TOl We find that there are two values of eV at which 
AIpAT changes drastically. First, there is a step located 
at eV = ££) and second, a Fano-like peak appears at 
eV — —€u- The (negative) derivative at eV = —cd 
seems to diverge. The high of both the peak and the 
step increase in a monotonous manner as a function of 
the ratio of the device voltage eVd divided by the dot level 
e D ■ When eVd is small, the peak is much higher than the 
step. Increasing eVd, the peak further increases, but the 
step height increases faster, starting from the threshold 
device voltage eVd = A — e^j. In Fig. [TOl we find that for 
f-D ~ 0.4 and eVd = 0.8, the peak height is still higher 
than the step, but with = 0.6, and eVd = 0.6, the 
step becomes higher than the peak. Here for specificity, 
we only consider the case where eo > 0, but results for 
££) < can be obtained in a similar manner, exploiting 
electron hole symmetry. 

In order to further understand the behavior of AIpat, 
we consider the different contributions of this photo- 
assisted current, which are shown in Fig. 1111 Specif- 
ically, we plot the elastic current rcnormalizcd by the 
environment, as well as the right and left inelastic cur- 
rents. We find that the elastic part is symmetric between 
eV positive and negative. It is almost equal to zero when 
\eV\ < CD- It shows a step at \eV\ = ed- This can be 
understood from the fact that at the threshold eV = eo 
electrons tunnel from the normal metal lead to the super- 
conductor predominantly by making resonant transitions 
through the dot. Electrons easily tunnel to the quantum 
dot, in a sequential manner becoming a Cooper pair in 
the superconductor. For eV = —eo, the same reasoning 
can be made for incoming holes, or equivalently for elec- 
trons exiting the superconductor: a Cooper pair in the 
superconductor is split into two electrons, which tunnel 
to the quantum dot and then to the normal metal lead. 
Because of the energy conservation condition, it is then 
necessary to have eV < —eo- 
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Turning now to the inelastic current, we consider the 
contribution of 2 electrons being transfered inelastically 
from the normal metal lead to the superconductor, which 
is called the J™"^' in Fig. [TTJ When eV < eo, electrons 
tunnel through the quantum dot to the superconductor 
by absorbing or emitting energy to environment. How- 
ever, as in the elastic case, the transfer from the normal 
metal to the dot is more favorable when eV > eo, which 
explains the presence of the step in at this bias volt- 



age. 

Next, we consider the contribution of two electrons 
tunneling from the superconductor to the normal metal 
lead, which is called J™'^' in Fig. [TlJ When eV is pos- 
itive, /™^' is nonzero only when eVd > 2eV, this case 
corresponds to the process of the Cooper pair absorb- 
ing energy from environment to tunnel to the normal 
metal lead. If 2eo < eV^, a small step occurs (not 
shown) at eV = eo corresponding to the activation of 
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FIG. 10: AIpAT plotted as a function of the detector bias voltage, for dot energy level: ed = 0.4 (left panel), 0.6 (right panel) 
and for several values of device bias eVd- 
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FIG. 11: Elastic and inelastic contributions to AIpat, as a function of detector bias voltage. The dot energy level is fixed at 
eo = 0.6. The mesoscopic device bias eVd is: 0.8 (left panel), 0.4 (center panel), and 0.3 (right panel). 



the two Cooper pair electrons on the dot: this small fea- 
ture can only be seen by zooming in the picture. When 
eV is negative, the absorption of energy from the envi- 
ronment becomes much more favorable. Because of this, 
the analog of the step corresponding to eV = — e_D hi the 
elastic current is smoothed out, and it saturates around 
eV = — e_D. Nevertheless, J™"^' also contains contribu- 
tions where electrons emit energy to the environment. 
Starting from eV < 0, the process of the Cooper pair 



tunneling to the normal metal lead and emitting energy 
to environment has first a small contribution to /™^' but 
it really becomes noticeable below eV = —sd and eventu- 
ally saturates for lower bias voltage (not shown) . As the 
voltage of the mesoscopic device is lowered (left to right 
panels of Fig. [TT|) . two things occur: first, the ampli- 
tude of all current contributions decreases, second, the 
smoothing of the J™"^' is reduced because the range of 
energy available to absorption and emission is reduced. 
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In brief the sum of the contributions for emission and 
absorption in /™^' have a tendency to compensate the 
elastic current at voltages where saturation is reached. 
^From the above considerations, we therefore can inter- 
pret the curve of IS.Ipat, and we understand when the 
detector circuit absorbs or emits energy from/to meso- 
scopic device: when eV is negative, the PAT current is 
mainly due to absorption for \eV\ < eu, and emission for 
\eV\ > eo- 

Note that in all our numerical results, the PAT current 
is in units of 2e^R'^jf-/^T{l - T)/Tr^h'^ A^Rk ■ To check 
the observability of these predictions we estimate: 71 = 
0.2A, 72 = 0.02A^ = 240/iel/, R = OMRk, T = 0.5 
(see also Ref. [Mnilillll) • This implies, e.g. for the 
PAT current in Fig. [10] when the detector bias is close 
to ±e_D, IpAT — 5 — lOpA with the device bias Vd = 
0.8A/e — iSfiV. This value is acceptable if one compares 
it with the value of current which has been estimated in 
Ref. [13. It is also acceptable with present day detection 
techniques. 



IV. CONCLUSION 

In conclusion, we have presented a new capacitive cou- 
pling scheme to study the high frequency spectral density 
of noise of a mesoscopic device. As in the initial pro- 
posal of Ref. the effect of the noise originating from 
a mesoscopic device triggers an inelastic DC current in 
the detector circuit. This inelastic contribution can be 
thought as a dynamical Coulomb blockade effect where 
the phase fluctuations at a specific junction in the detec- 
tor circuit junction are related to voltage fluctuations in 
same junction. In turn, such voltage fluctuations orig- 
inate from the current fluctuations in the nearby meso- 
scopic device, and both are related by a trans-impedance. 
The novelty is that here, because the junction contains 
a superconducting element, in the subgap regime, two 
electrons need to be transferred as the elementary charge 
tunneling process. In a conventional elastic tunneling sit- 
uation, the two electrons injected from (ejected in) the 
normal metal lead need to have exactly opposite energies 
in order to combine as a Cooper pair in the supercon- 
ductor. Here, this energy conservation can be violated in 
a controlled manner because a photon originating from 
the mesoscopic circuit can be provided to/from the con- 
stituent electrons of the Cooper pair in the tunneling 
process. 

We have computed the DC current in the detector cir- 
cuit for two different situations. In a first step, we con- 
sidered a single NS junction, and we computed all lowest 
order inelastic charge transfer processes which can be in- 
volved in the measurement of noise: the photo-assisted 
transfer of single (and pairs of) electrons (with energies 
within the gap) into quasiparticle(s) above the gap, and 
photo-assisted Andreev transfer of electrons as a Cooper 
pair in the subgap regime. It was shown that the latter 
process dominates when the source drain bias is kept well 



within the gap. Close but below the gap, the absorption 
of quasiparticle dominates, and we observe a crossover in 
the current between the Andreev and quasiparticle con- 
tributions. For the above processes, we demonstrated 
the dependence of the DC current on the voltage bias 
of the mesoscopic device, chosen here to be a quantum 
point contact. When this bias eVd is increased, the over- 
all amplitude of the spectral density and its width scale 
as Vd, so that the phase space (the energy range) of elec- 
trons which can contribute to the Andreev processes is 
magnified. We therefore have gained an understanding 
about how the measurement of the DC current can pro- 
vide useful information on the noise of the mesoscopic 
device. Nevertheless, one should point out that with this 
NS setup, it is difficult to isolate the contribution of the 
photo-assisted current which involves, respectively, the 
absorption and the emission of photons from the meso- 
scopic device. For biases close to the chemical potential 
of the superconductor, both will typically contribute to 
the photo-assisted current. 

Next, we considered a more complex detector circuit 
where the normal metal and the superconductor are sep- 
arated by a quantum dot which can only accommodate 
a single electron at a given time. There, the dot acts 
as an additional energy filtering device, with the aim of 
achieving a selection between photon emission and ab- 
sorption processes. We decided to restrict ourselves to 
the photo-assisted Andreev (subgap) regime, assuming 
that the dot level is well within the gap. By comput- 
ing the total excess photo-assisted current as well as its 
different contributions for absorption and emission, and 
right and left currents, we found that for DC bias volt- 
ages comparable to eV = —cd it is possible to make such 
a distinction. The NDS detection setup could therefore 
provide more information on the spectral density of noise 
than the NS setup, but its diagnosis would involve the 
measurement of smaller currents than the NS setup, be- 
cause of the presence of two tunnel barriers instead of 
one. 

Note that normal lead - quantum dot - superconduc- 
tor setups have already been investigated theoreticalljs^^ 
and experimentally recently^^. In such works, the em- 
phasis was to study how the physics of Andreev refiec- 
tion affected the Kondo anomaly in the current voltage 
characteristics. In Ref. [2^, the quantum dot consisted of 
a carbon nanotube making the junction between a nor- 
mal metal lead and a superconductor. Here, we did not 
consider the quantum dot in the Kondo regime, and we 
included interactions on the dot in the Coulomb blockade 
regime. 

A central point of this study is the fact that all con- 
tributions to the photo- assisted current, for both the NS 
and NDS setups, can be cast in the same form: 

AIpATieV) « jdn^-^^^StA±^)Kprocess{^,eV, ...)(52) 

Where Kprocess is a Kernel which depends on the na- 
ture (elastic or inelastic) as well as the mechanism (sin- 
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gle quasi-particle or pair tunneling) of the charge trans- 
fer process. When deaUng with an elastic process, one 
understand that the environment renormalizes the DC 
current even when no photon are exchanged between the 
two circuits. In the case of inelastic tunneling only, the 
frequency Q corresponds to the total energy of the two 
electrons which enter (exit) the superconductor from (to) 
the normal metal lead. Finally, the sign of the frequency 
(and thus the bound of the integral in Eq. ((52)) . which 
are left "blank" here) decides whether a given contribu- 
tion corresponds to the absorption or to the emission of 
a photon from the mesoscopic circuit. 

The present proposal has been tested using a quantum 
point contact as a noise source, because the spectral den- 
sity of excess noise is well characterized and because it 
has a simple form. It would be useful to test the present 
model to situations where the noise spectrum exhibits 
cusps or singularities. Cusps are known to occur in the 
high frequency (close to the gap) noise of normal super- 
conducting junctions. Singularities in the noise are know 
to occur in chiral Luttinger liquid, tested in the context 
of the fractional quantum Hall effect^. Such singulari- 
ties or cusps should be easy to recognize in our proposed 
measurement of the photo-assisted current. 

On general grounds, we have proposed a new mech- 
anism which couples a normal metal - superconductor 
circuit to a mesoscopic device with the goal to under- 
stand the noise of the latter. The present setup suggest 
that it is plausible to extract information about high fre- 
quency noise. High frequency noise detection now consti- 
tutes an important subfield of nano/mesoscopic physics. 
New measurement setup schemes which can be placed 
"on-chip" next to the circuit to be measured are useful 
for further understanding of high frequency current mo- 
ments. 
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V. APPENDIX A 



In this Appendix, we define 4*0^, i^|e^(w, eV, 77), and 
K^'it{^,eV,'n) in the 



-A-eV 



lA+eV 

+ IT] \d + irj — trj J 



A-eV 
A+eV 

1 1 \ / 1 1 



X 2- 



5 + ir] S + uj + ir] J \S + it] S — iij 



S + 17] \6 — Lu + irj 6 + Lu — irj 



pQ.~2A-2eV 

KrMn,eV,7j) = / d6 

J2A+2eV-0 



(54) 



X ^/{n + S- 2ey)2 - 4AV(f^ -S- 2eVy - 4A2 
/ 1 1 \ 

X 



(55) 



VI. APPENDIX B 

In this Appendix, we recall the definition of Keldysh 
Green's functions^l. First, we define the anomalous 
Green's function, describing the pairing of electrons with 
opposite spin in the superconductor: 

F,{q,t-t') = -{TKC-q,-a{t)Cq,a{t')) , 



{TKclAt)cl 



r{t')) 



If both t and t' are in the upper branch, and t > t' 
or both t and t' are in the lower branch, and t' > t 
then F,{t+-t'+) = Fiit^-t'_) ^ sgnia)ugVge-'^^^'~''\ 
These Green's functions enter the description of the An- 
dreev process. 

If we consider the single quasiparticle tunneling in the 
superconductor, we use the conventional definition of the 
Green's function as for normal metals. 

Secondly, we define the Green's function of the one 
level QD: Goait - t') = {TKCDa{t)c+^{t')) . 

Simplifications occur because the quantum dot has a 
singly occupied level with energy e/j, the first electron 
is transferred to the lead before the second hops on the 
quantum dot so that in our work, we only consider the 
QD Green's function where both time quantities t and t' 
arc in the upper or lower branch, and the Green's func- 
tion values only when t > t' ii t, t' in the upper branch 
G^^(i - t') = e"'"" and t' > t \i t, t' in the lower 
branch, then G*jy^{t - t') = e-"o(*-*'). 

The Green's function in the normal metal lead reads: 

GLa{k,t-t')={TKCka{t)cl„{t')). 

In our work, we consider the cases where two electrons 
tunneling from/to normal metal lead, so that we only 
consider normal metal Green's function where t and t' are 
in the different branches. For the case electrons tunneling 
from the superconductor to the normal metal, we use the 
greater Green's function G^^{k,t-t') = e-^('='=-''i)(*-*') 
with Efc > /ii. For the case electrons tunneling from the 
normal metal to the superconductor, we use the lesser 
Green's function G^^(fc,t-t') = -{cl^{t')cka{t)) = 
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VII. APPENDIX C 

111 this part, we present the denominator products 
which appear in the tunnehng current through the NS 
junction as a Cooper pair. Such denominators come from 
the energy denominators of the transition operator T. 
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(56) 
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Change variables in I™*^' , /™^' as 



5 = e-e' . 



Wc define 



Ylx.ri) = / dE , — 



(57) 



(58) 



7r + 2arcsin(^) 
2yA2~~(^TI^ ' 



then we define the weight functions as 



eV 

de { [2x(-e, -rj) + x{-f- - ^, -v)] [x(e, '7) + x(-e, '?)] 



+ x(-e, ~V) [x(e - 77) + x(-e - t^, v)]} 



(60) 



if™^'(r!, 6^,77) = / 

J2e 



f2-2eV 

3V-n 
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VIII. APPENDIX D 
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In this Appendix, we present the denominator prod- 
ucts which appear in the tunnehng current through the 
NDS junction. 

D'^ is the original denominator which is not affected 
by the environment 
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and is the denominator product affected by the environment 
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(63) 



where E)™^^ is the denominator product attributed to the inelastic current (affected by environment) and it is defined 
as 



1 



1 
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{eo - e' + «??)(e_D ~ e - iri) 
1 
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1 1 
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Since D°_^{e) = Dt{-e), Dt{e) = Dt{-e), so ^''^(e) = 
*K-e), *!i(e) = l'^(-e). But in fact, if we change the 
name of variable e —^ e' then change variable e' ~ — e, we 
will obtain the same formula for both cases eV > and 
eV < 0. Since e, e' arc independently equivalent, so that 
it is evidently ^'™'='(e,e', en, 77) = (e, e', e^,, 77). Here- 
after, we neglect the <— or ^ index in these functions. 

If eV > 0, the elastic current contributions in exist 
but the clastic current contributions in vanish (in 
contrast to the case of eV < 0). 

We change variables in inelastic contributions as 



and define 

K'r^Dsi^, eV, ep, 77) = /''de1'^'-*°*(e, e^,, 77) , (69) 
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